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We advocate charmed-hadron inclusive hadroproduction as a laboratory to probe in- 
trinsic charm (IC) inside the colliding hadrons. Working at next-to-leading order in 
the general-mass variable-flavor-number scheme endowed with non-perturbative frag- 
mentation functions recently extracted from a global fit to e + e _ annihilation data from 
KEKB, CESR, and LEP1, we first assess the sensitivity of Tevatron data of D° , D + , 
and D* + inclusive production to the IC parameterizations provided by Pumplin et al. 
We then argue that similar data from pp collisions at RHIC would have the potential to 
discriminate between different IC models provided they reach out to sufficiently large 
values of transverse momentum. 

1 Introduction 

The general-mass variable- flavor-number scheme (GM-VFNS) provides a rigorous theoretical 
framework for the theoretical description of the inclusive production of single heavy-flavored 
hadrons, combining the fixed-flavor-number scheme (FFNS) and zero-mass variablc-flavor- 
number scheme (ZM-FVNS), which are valid in complementary kinematic regions, in a uni- 
fied approach that enjoys the virtues of both schemes and, at the same time, is bare of their 
flaws. Specifically, it resums large logarithms by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution of non-perturbative fragmentation functions (FFs), guarantees 
the universality of the latter as in the ZM-VFNS, and simultaneously retains the mass- 
dependent terms of the FFNS without additional theoretical assumptions. It was elaborated 
at next-to-leading order (NLO) for e + e~ annihilation [2 photo- [3] and hadroproduction 
[H [5] . In this presentation [T] , we report recent progress in the theoretical description of 
£>-meson hadroproduction and argue that this process provides a useful probe of IC in the 
proton In Sec. [21 we present updated predictions for the transverse-momentum (pt) 
distributions of I? , D + , and D* + mesons [6] using our new FF sets [2]. In Scc.[3l we briefly 
review the constraints on IC obtained through an extension [7] of the CTEQ6.5 global anal- 
ysis [5] endowed with three representative IC models. In Sec. [H we investigate, for D° 
hadroproduction at the Tevatron and at RHIC, the shifts in cross section produced by these 
IC implementations. Our conclusions are summarized in Sec. [3] 



2 .D-meson hadroproduction revisited 

Our previous NLO predictions for inclusive D-meson hadroproduction [5] were evaluated 
with the FF sets of Ref. [5] . These and previous FF sets [TU] had been determined through fits 
to LEP1 data [11] in the ZM-VFNS. These determinations have recently been updated and 
improved [2] by also including KEKB and CESR HU data, switching to the GM-VFNS, and 
also incorporating hadron-mass corrections. They greatly improve the agreement with the 
Pt distributions of D°, D + , and D* + mesons measured at the Tevatron [13] , as may be seen 
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in Fig. [lla)-(c) and (d)-(f) on an absolute scale and in the data-over-theory representation, 
respectively. 
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Figure 1: The solid lines represent the default predictions, and the dashed lines indicate the 
theoretical uncertainty. 



3 Intrinsic charm from global analysis 

In Ref. [7j, the CTEQ6.5 global analysis [B] was repeated allowing for IC according to the 
Brodsky-Hoyer-Peterson-Sakai (BHPS) [13], meson-cloud [TS], and sea-like [7] models, with 
a modest and a marginal value of (x) c+ c in each case. The x distributions at fi = 2 GeV 
of the resulting c and c parton distribution functions (PDFs) are shown in Fig. HUa)-(c) for 
the BHPS, meson-cloud, and sea-like models, respectively. 

4 Probing intrinsic charm at Tevatron and RHIC 

In the GM-VFNS, D- meson hadroproduction receives sizeable contributions from partonic 
subprocesses with incoming c or c quarks [U [5] and is thus sensitive to their PDFs and 
IC- induced enhancements in them. In Fig. [3ja), the relative enhancements due to the six 
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Figure 2: The lower and upper dashed lines correspond to IC with modest and marginal 
(x) c+ c values, respectively, and the solid ones to zero IC. 



IC implementations of Ref. [7] are compared with the D° data of Ref . [13] . We observe that 
the two large-pr data points tend to disfavour the sea-like IC with marginal value of (x) c +c, 
although a firm statement would be premature. On the other hand, the effects due the IC 
of the BHPS or meson-cloud models are too feeble to be resolved by the presently published 
CDF data. 

Since IC typically receives a large fraction x of momentum from the parent hadron, 
especially in the BHPS and meson-cloud models, and the kinematical upper bound on x 
scales with xt = %Pt/\/~s, the sensitivity to IC may be enhanced by increasing pr or 
lowering the cm. energy yfs. Keeping in mind that the CDF analysis of Ref. [13] is merely 
based on 5.8 pb _1 of data recorded in February and March 2002 and that, since then, the 
integrated luminosity of run II has increased by more than a factor of 1000, exceeding 6 fb _1 
as it does, significantly more precise data at larger pr values are expected to come, with an 
increased sensitivity to IC, as may be seen from Fig. EJb). A hadron collider with smaller 
cm. energy is in operation, namely RHIC with presently y/s = 200 GeV, where the px 
distribution of pp — » X c + X could be measured. So far, only the STAR Collaboration has 
published Z3-meson production data, namely for dAu — ► D° +X with pt < 2.5 GeV [16]. In 
Fig. [3Jc), we observe significant enhancements due to IC in pp collisions at yfs = 200 GeV, 
by up to 75% at p T = 10 GeV. The high-energy pp mode of RHIC, with y/s = 500 GeV, will 
accrue more luminosity, perhaps 200-400 pb" 1 . However, this will happen at the expense 
of lowering the accessible x values and, thus, of reducing the relative IC effects on the cross 
sections, especially for the BHPS and meson-cloud models, as is evident from Fig.[3jd). 

5 Conclusions 

In the case of the Tevatron, we found that the BHPS and meson-cloud models yield insignif- 
icant deviations from the zero-IC predictions, while the shift produced by the sea-like model 
is comparable to the experimental error and tends to worsen the agreement between theory 
and experiment. However, the experimental errors in Ref. |13] are still too sizeable to rule 
out this model as implemented in Ref. [7] . This is likely to change once the full data sample 
of run II is exploited. 

As for pp collisions with yfs = 200 GeV at RHIC, all three IC models predict sizeable 
enhancements, by up to 75% at px = 10 GeV. In a future high-energy pp-collision mode of 
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RHIC, with -y/s = 500 GeV, which is to accrue more luminosity, the cross sections would be 
increased, while the relative shifts due to IC would be considerably smaller for the BHPS 
and meson-cloud models. 
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Figure 3: The solid (dashed), densely dotted (dot-dashed), and scarcely dotted (dotted) 
lines correspond to the BHPS [14], meson-cloud [15], and sea-like [7] models with modest 
(marginal) values of {x) c +c- 
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